Resonant electronic states and I-V curves of Fe/MgO/Fe(100) tunnel junctions 
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The bias dependence of the tunnel magnetoresistance (TMR) of Fe/MgO/Fe tunnel junctions 
is investigated theoretically with a fully self-consistent scheme that combines the non-equilibrium 
Green's functions method with density functional theory. At voltages smaller than 20 mVolt the 
I-V characteristics and the TMR are dominated by resonant transport through narrow interface 
states in the minority spin-band. In the parallel configuration this contribution is quenched by a 
voltage comparable to the energy width of the interface state, whereas it persists at all voltages in 
the anti-parallel configuration. At higher bias the transport is mainly determined by the relative 
positions of the Ai band-edges in the two Fe electrodes, which causes a decrease of the TMR. 
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Modern magnetic sensors, such as read heads for hard 
disk drives, are based on the tunnel magnetoresistance 
(TMR) effect. This is the drop in resistance of a magnetic 
tunnel junction (MTJ) formed by two magnetic layers 
when the mutual alignment of their magnetization vec- 
tors changes from antiparallel (AP) to parallel (P). The 
TMR magnitude is given by TMR = [J P - Zap] /.Zap, 
with I a (V) being the current at the voltage V for the 
a configuration (P or AP). Huge TMR ratios have been 
recently achieved in epitaxial, all crystalline Fe/MgO [l| 
and CoFeB/MgO 0] MTJs, reaching up to 500% at room 
temperature and 1010% at 5 K 0. These large values of 
TMR are broadly attributed to the phase coherent and 
transverse momentum conserving transport. The TMR 
thus is governed not only by the spin-polarization of the 
electrode density of states (DOS), but also by the details 
of the wave-functions matching across the barrier. Mag- 
netic transition metals with bcc crystal structure (Fe) 
possess a high transmission A i band, which decays slowly 
across the MgO barrier This is fractionally filled 
for majority spins (f), and empty for minority spins (J.). 
Since these bands dominate the tunneling current, bcc 
transition metals with MgO barrier effectively behave as 
half-metals, and the TMR is expected to be very large 

an. 

Another important, but much less investigated aspect, 
is the relation between electronic states and the I-V 
characteristics in these highly crystalline MTJs. Inter- 
facial states and details of the Fe band-structure, oth- 
erwise washed out by disorder, play an important role 
in the transport and indeed can be identified through 
the I-V curves and its derivatives (G(V) = dl/dV and 
S(V) = d 2 I/dV 2 ). For instance, high-quality conven- 
tional MTJs (2-3 nm MgO thickness) show a pronounced 
broad peak on the S(V) curve at about 1 Volt for the 
AP, and a number of small peaks at lower voltages in 
the P configuration 0, 0] ■ Combined with a quantitative 
theory these measurements can provide a wealth of infor- 
mation, and help on the level of device design. We note 



that very thin MTJs, with a MgO thickness of about 
1 nm, have to be used for ultra high recording density 
(>500 Gbit/inch 2 ) HDD readers [8J. Remarkably, at 
these small MgO thicknesses, the growth mode of the 
Fe/MgO/Fe stack changes 0- 

In this letter we investigate theoretically the I-V char- 
acteristics of perfectly crystalline Fe/MgO MTJs, and 
demonstrate that its features originate from the sweep- 
ing of the Ai band-edges and of interface states across 
the bias window. These generic features are emphasized 
here, as they resemble closely those studied in molec- 
ular devices Q and magnetic point contacts (ic| . We 
investigate the ultra-thin MTJ regime, with the goal 
to provide a solid basis for decoding future I-V mea- 
surements of < 1 nm thick MgO barriers, in term of 
their relation to the underlying electronic structure and 
possibly growth defects. Calculations are performed us- 
ing SMEAGOL [11], which combines the non-equilibrium 
Green's function method with density functional theory 
(DFT) [13]. The total transmission coefficient T(E; V) 
is self-consistently evaluated at finite bias and integrated 
to give the spin-current. 



(1) 



Here a is the spin index (T, I) and f± is the Fermi func- 
tion calculated at (E — Ep ± eV/2)/kBT, with E-p the 
Fe Fermi energy, &b the Boltzmann constant and r the 
electronic temperature. Translational invariance allows 
to write T a {E; V) = ± J T a (E, kj_; V) dk ± , where the 
integral runs over the 2D Brillouin zone perpendicular 
to the transport direction with area £lk- Large k^-point 
samplings are necessary to converge T(E\ V), so that an 
extremely stable numerical algorithm is needed [HI]. 

We consider a MTJ oriented along the Fe(100) direc- 
tion and formed by 4 MgO monolayers (~1.1 nm long) . 
The atomic coordinates are those from reference [14] . 
which have been optimized by high-accuracy plane-wave 
calculations. The unit cell used for the transport con- 
tains also 8 Fe atomic layers on each side of the MgO 
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FIG. 1: (Color online) (a) Planar average AVh of the differ- 
ence between the Hartree potential at 0.5 Volt and the one 
at 0-bias. (b) Planar average Ap of the difference between 
the charge density at 0.5 Volt and the one at 0-bias. The 
diamonds and dots indicate the location of the Fe and MgO 
layers. 



barrier, which are enough to screen the charge density 
from the interface. A 7x7 k^-points mesh converges the 
charge density, but a 100x100 mesh is used for evaluat- 
ing T(E;V). The basis set employed is single-^ for the 
Fe p and d orbitals, while double-C is used for all other 
orbitals [l2j]. The local density approximation (LDA) is 
adopted throughout, and we use a real space mesh cutoff 
of 600 Ry and r = 300 K. 

First, we look at the electrostatic Hartree potential 
drop across the junction. In Fig. [Ha) we show the 
difference AVh between the planar average of the self- 
consistent Hartree potential at a finite bias and that at 
0-bias along the junction stack (z-axis). AVh is flat in the 
electrodes and decays linearly in the MgO, demonstrating 
that the usual approximation of shifting rigidly the elec- 
tronic structure of the electrodes and then applying a lin- 
ear potential drop across the barrier [la ] is well justified. 
In Fig. [Tfb) we also show the difference Ap between the 
planar average of the self-consistent charge density at a fi- 
nite bias and at 0-bias, so that Ap(z) oc — d 2 AVn(z) / dz 2 . 
At all voltages we find Ap(z) increasing linearly with V, 
and charge accumulating at the extremal layers of the 
electrodes just before MgO. Inside the MgO Ap(z) oscil- 
lates due to the electric field induced polarization. This 
is confirmed by a DFT calculation for an isolated MgO 
slab of the same thickness in an equivalent electric field, 
which shows analogous oscillations. 

We start our analysis of the transport properties by 
presenting T(E; V) (Fig. [2]). Three main features appear 
in the V = transmission coefficient: © a sharp increase 
(note the logarithmic scale) in transmission at around - 
1 eV for the j spins in the P configuration, © a similar, 
although smoother increase above +1 eV for the J. spin 
in the P configuration and for the AP configuration, and 
® a sharp resonance at Ep for the j spins in the P con- 
figuration, which is still visible in the AP configuration. 

The first two features are associated with the Ai band- 
edges, respectively for the f and { spins, as it can be 
clearly seen in Fig. [3J Note also that these band-edges 
coincide with energy regions where the average number 
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FIG. 2: (Color online) Transmission coefficient T(E; V) as 
a function of energy, E, and for different biases, V. The 
vertical lines are placed at E = Ef ± eV/2 and enclose the 
bias window. The dashed (solid) line is for the j (J.) spin 
band in the P configuration, while the dashed-dotted line is 
for the AP configuration. The lowest panel is T(E; 0) for the 
AP configuration and different broadening S. 
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FIG. 3: (Color online) (a) Bulk Fe band-structure along the 
r — > H direction (the bands with Ai-symmetry are empha- 
sized), (b) T(E;0) in the P configuration, (c) average number 
of channels per fc-point for bulk Fe, n c , (d) bulk Fe DOS, (e) 
interface Fe-layer DOS. Note that the Ai band-edges coincide 
with a sharp increase in the transmission coefficient T(E; 0) 
and with a large n c . 



of Fe-channels per k^-point n c is maximized [Fig. G2c)]. 
Since these are a feature of the electronic structure of 
Fe alone, their position as a function of bias is set by 
the quasi Fermi energies, E-p ± eV/2, of the two elec- 
trodes. For instance the sharp increase of T^(E;V) in 
the P configuration (©) moves to -0.75 eV and -0.25 eV 
respectively for voltages of 0.5 Volt and 1.5 Volt, follow- 
ing the a\ band-edge of the left lead. In the same way 
the broad peak in both P and AP at about 1.5 eV above 
Ep (@) splits into two broad sub-peaks, separated by 



FIG. 4: (Color online) Schematic representation of the bias 
dependence of surface states localized at the interfaces be- 
tween the magnetic electrodes and the tunnel barrier, (a) 
and (b) are for the left (L) and right (R) interface states hav- 
ing the same energy and therefore resonating at V=0; (c) and 
(d) are for the L and R interface states not resonating at 
zero-bias. 
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FIG. 5: (Color online) Current I, S(V) and TMR as function 
of voltage V for different values of the imaginary part of the 
energy S: (a) I-V for P, (b) I-V for AP, (c) S{V) for P, (d) 
S(V) for AP, (e) TMR. In the inset of panel (e) we show the 
TMR in the low-bias region for 5 = 0. 



eV. Note however that whereas the height of the peak is 
roughly constant for the AP configuration, in the P con- 
figuration the sub-peak entering the bias window shrinks 
drastically. 

In contrast the peak in transmission at Ep for the J, 
spins in the P configuration and in the AP configuration 
cannot be associated with leads band-edges but instead 
originates from a narrow interface state in the minority 
band 

[IB El- This 

is spatially localized at the inter- 
face between Fe and MgO and it is resonating at Ep. A 
comparison between the bulk Fe DOS and the DOS for 
the interface Fe layer [Figs. [Sfd) and (e)] indicates the 
surface state for J, spins at Ep [4| . The dynamics of such 
a interface state under bias can be understood from the 
cartoon of figured] Consider the panels (a) and (b) where 
two identical interface states are localized on either sides 
of the tunnel barrier. This is the situation encountered 
here for the J, spins. The transport is then resonating 
across the barrier at the interface state energy Er = E^. 
In the case considered here the resonance energy is close 
to E F . 

Since, in general the state is coupled more strongly to 
one of the electrodes, it will trace closely its quasi Fermi 
energy. For instance for positive bias and a interface 
state localized on the left-hand (right-hand) side of the 
junction [Fig. S[b)] we obtain E L (V) = E L (0) + eV/2 
(Er(V) = £r(0) - eV/2). This brings the states on 
either side of the junction out of resonance and generally 
suppresses the transmission. Thus the peak in T(E;0), 
originating from a resonating interface state across the 
barrier, will evolve into two smaller peaks separated by 
an energy eV. Indeed this is the behavior observed in 
Fig. O For instance the two peaks centered at Ep are 



separated by 0.5 eV and 1.5 eV respectively for voltages 
of 0.5 Volt and 1.5 Volt. In the AP configuration this 
situation is found at any bias, since the spin J. interface 
state is always present on one side of the junction only. 

A second possible situation is when the interface states 
on the left-hand and right-hand side of the tunnel junc- 
tion have a different origin and are placed at different en- 
ergy. In this case we do not expect zero-bias resonance, 
however there will be a critical voltage at which the res- 
onant condition is met. In this case we expect the rise 
of a large peak in the transmission coefficient at a bias 
eV = Ei,(0) — £Jr(0). This situation has never been en- 
countered for the symmetric MTJ investigated here, but 
it is likely to be the most typical case in real junctions. 

We now move to analyzing the I-V characteristic and 
the TMR (Fig. O, starting from the low bias region (V < 
0.4 Volt). The most apparent feature in this bias range 
is a sharp reduction of the TMR from its zero-bias value 
followed by a rapid increase which peaks at V » 0.3 Volt. 
The sharp reduction can be associated with a decrease in 
the P I current originating from the loss of the inter- 
face state resonant condition at 7 « 20 mVolt, a bias 
which roughly corresponds to the line-width of the inter- 
face state. For V < 20 mVolt the P current is shared 
by the two spin-species, while for V > 20 mVolt the f 
component dominates. Such a reduction in the P j com- 
ponent at the bias corresponding to the resonant condi- 
tion loss can be clearly observed in the S(V) plot of Fig. 
E£c). As the bias further increases the conductance in 
the P configuration is approximately constant. In con- 
trast the conductance of the AP configuration is reduced 
for V between 0.1 and 0.3 Volt, and this behavior results 
in an increase of the TMR with the broad peak at 0.3 
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Volt. 

The high bias region is characterized by a decrease 
of the conductance in the P configuration at around 
1 Volt and a dramatic increase of the AP current for 

V > 1.5 Volt. This produces a strong reduction of the 
TMR with bias and an almost complete suppression for 

V > 1.5 Volt. In fact at about 1.75 Volt the TMR be- 
comes negative. Such a voltage range should be put in 
comparison with the band offset. In our LDA calcula- 
tions E-p for Fe is positioned ~1.8 eV below the MgO 
conduction band minimum and ^3.0 eV above the va- 
lence band maximum. This means that voltages of the 
order of 1.5 Volt are still rather far from those needed for 
tunneling across a reduced barrier. 

The high-bias behavior is dominated by the relative 
energy shift with bias of the Fe states. The origin of 
the reduction of the conductance in the P configuration 
for V > 1 Volt is that once the A{ band-edge of the left 
lead enters the bias window, the high transmission re- 
gion only extends over part of such window (Fig. [2]), so 
that the increase of the current with bias is reduced by a 
factor of about 2. Note that the current in the J, is neg- 
ligible. The broad peak of the J, transmission at about 
1.5 eV above the Fermi energy, related to the A} band- 
edge, never contributes significantly to the current in the 
P configuration. The reason is that whereas for positive 
bias the A{ band-edge of the right lead enters the bias 
window, the one on the left lead moves away from the bias 
window, so that the peak shrinks as bias is applied. The 
general evolution with bias of this transmission peak re- 
sembles the one of the surface states schematically shown 
in Figs. HUa) and (b). In the AP configuration however 
this peak does contribute to the current, since electrons 
belonging to the A{ band in the right lead that are inside 
the bias window can tunnel into the A{ bands of the left 
lead. This leads to an increase of the conductance once 
the Aj band-edge moves into the bias window, resulting 
in a drastic increase of the AP current at V > 1.5 Volt. 
We note that our results in the high bias region are in 
contrast with those in reference [18j |. in which a rapid 
increase of the conductance is found also for the P con- 
figuration. 

It is also interesting to comment on the S(V) plots 
[FigEJc) and (d)]. For both the magnetic configurations 
one can observe a peaked structure. This is observed also 
experimentally Q and attributed to resonances in the 
transmission. At the qualitative level this interpretation 
is confirmed by our results in the low bias region, where 
the interface states determine the behavior of the I-V. 
The actual position of the peaks however depends on 
the atomic details of the two interfaces and we do not 
expect agreement. The high-bias region however is more 
problematic. In particular we note that our I-V for the P 
configuration does not increase as rapidly as that found 
in typical experiments. One possible explanation for such 
difference is that the MgO barrier calculated here is very 



thin and the current consequently already large at small 
bias, another reason could be presence of defects at the 
interface and in the MgO in experiments. 

Finally we investigate the effect of generic disorder, 
motivated by observations 0] • This is modeled at a sim- 
ple level by adding a small imaginary part 6 to the energy 
when calculating T(E; V), i.e. it corresponds to a uni- 
form level broadening [17| . The effects of such addition 
on the I-V, S(V) and TMR are shown in Fig. El while 
those on T(E; 0) for the AP configuration are shown in 
Fig. O Fig. [5] clearly shows that increasing S results in a 
gradual suppression of the T(E; 0) resonance at E-p. For 
6 = 10~ 4 Ry the peak is still clearly visible, and only for 
8 = 1CP 3 Ry it completely disappears. As a result the 
TMR at low bias largely increases and the non-monotonic 
behavior for V < 0.4 Volt is suppressed. In contrast the 
high bias region is barely affected by S. These results, al- 
though indicative of the effects of disorder on the TMR, 
should be taken with caution. The broadening S intro- 
duces unstructured disorder, and the transmission of all 
the spin-channels is equally reduced. In reality one may 
expect the transmission to either increase or decrease de- 
pending on the type of scattering center, which in general 
will act differently on the different spin-channels. This 
will in general result in an enhancement of the current in 
the AP alignment, causing a reduction of the TMR. 

In conclusion, we investigated the bias dependence of 
the TMR of an epitaxial Fe/MgO/Fe tunnel junction. 
We identify two different bias regions, which are affected 
by two types of electronic resonances: (i) interface reso- 
nance states and (ii) band edges. The low bias region is 
characterized by a non-monotonic behavior of the TMR 
caused by resonant scattering across the barrier due to 
the highly localized interface states. In the high bias re- 
gion (V > 0.4 Volt) the TMR decreases monotonically, 
mainly due to band edge states. We have also shown 
how disorder can suppress the transport through inter- 
face states and smoothens the TMR. 
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